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ABSTRACT 


The longitude -dependent part of the geopotential usually gives rise only to short- 
period effects in the motion of an artificial satellite. However, when the motion of 
the satellite is commensurable with that of the earth, the path of the satellite repeats 
itself relative to the earth and perturbations build up at each passage of the satellite 
in the same spot, so that there can be important long-period effects. 

In order to take these effects into account in deriving a theoretical solution to the 
equations of motion of an artificial satellite, it is necessary to select terms in the 
longitude -dependent part of the geopotential that will contribute significantly to the 
perturbations. We have tried to make a selection that is valid in a general case, 
regardless of the initial eccentricity of the orbit and of the order of the resonance. 

The solution to the equations of motion of an artificial satellite, in a geopotential 
thus determined, is then derived by using Hori's method by Lie series, which, by its 
properties regar ding canonical invariance, has proved advantageous in the classical 
theory. 


RESUME 

La partie du geopotentiel , dependant de la longitude, ne provoque 
generalement que des effets a court e periode dans le mouvement d'un satellite 
ifici el . Cependant , quand le mouvement du satellite est commensurable 
avec celui de la terre, le parcours du satellite se repete par rapport a la 
terre et les perturbations s'ajoutent a chaque passage du satellite au meme 
endroit, et peuvent ainsi creer d'importants effets a longue periode. 
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Pour tenir compte de ces effets en derivant une solution theorique des 
equations du mouvement d'un satellite artificiel, il est necessaire de choisir 
ceux des termes de la partie du geopotentiel dependant de la longitude, qui 
contribueront d’une fagon significative aux perturbations. On a essaye de 
faire un choix qui soit valable pour le cas general, quelques soient la pre- 
miere excentricite de l'orbite et l'ordre et l'ordre de la resonance. 

On a ensuite derive la solution des equations du mouvement d'un satellite 
artificiel dans tin geopotentiel ainsi determine, en employant la methode de 
Hori par la serie de Lie, qui a ete prouvee avantageuse dans la theorie 
classique, grace a ses proprietes concernant 1' invariance canonique. 

KOHCITEKT 

3aBHCHina h ot flOJiroTbi cocTaBJunoiuan reonoTeHiiHaiia oCbWHO 
BbI3bIBaeT JIHUIb KOpOTKOnepHOUHbie H3MeHeHHH B flBHDKeHHH HCKyCCTBeH- 
hoto cnyTHHKa aeMJiH. OnHaKO, ecjiH flBMsceHHe h CK yccTBGHHoro 
enyTHHKa COH3MepHMO C HBHJKeHHeM 3 eMJTH , TO nyTL HCKy CCTBeHHOrO 
cnyTHHKa ot hochtcjib ho 3eMJiH h nepTypgauHfi, nocTpoeHHbix Ha 
KaiKflOM oTpe3Ke nyTH HCKy ccTBeHHoro cnyTHHKa b toh see Tonxe, 
noBTopaeTcn, TaK hto boshohctbhh yace MoryT Ghtb nojironepnon- 

HblMH . 


fljiH yneTa sthx bo3h9hctbhh npn nonyneHHH TeopoTHnecKoro 
pemeHHH ypaBHeHHH ubhscghhh HCKy ccTBeHHoro cnyTHHKa HeoCxoflHMo 
BbiCpaTB T© HJieHbi 3aBHCHiaeH ot nojiroTbi cocTaBJifliomeH reonoTeH- 
UHajia’, KOTopbie abjihiotch Ba>KHOH nacTbio nepTyp6au,HH. Mbi nbiTajincb 
ocymecTBHTb BbiCop TaKHM oC : pa30M, htoCh oh ro,i(HJicH mm oCidero 
cjiynan, b He3aBH chmocth ot HanajiBHoro OKcqeHTpHCHTeTa h nopajiKa 
pe30HaHca. 

3aTeM b reonoTeHiiHane, onpenejreHHOM TaKHM cnocoCoM, bm- 
bouhtch peine Hne ypaBHeHHH' nBHKeHHfl HCKy ccTBeHHoro cnyTHHKa no 
MeTOfly XopH c noMoiiibio cepnn JIh. 3tot MeTOfl 0 Ka 3 ajtcH noxesHbiM 
b KJiaccHnecKOH TeopHH CjiaronapH CBoefi ycTaHOBHBineHCH HeH3MeH— 


HOCTH. 



ON THE TESSERAL-HARMONICS RESONANCE 
PROBLEM IN ARTIFICIAL-SATELLITE THEORY 

Barbara A. Romanowicz 
I. INTRODUCTION 

If the gravitational potential of the earth is expanded in terms of Legendre poly- 
nomials and functions, then in order to obtain a good approximation in the determination 
of the orbit of an artificial satellite, it is usually sufficient to consider the zonal, lon- 
gitude-free terms of the expansion. However, the influence of the tesseral terms 
becomes important in the case when the mean motion of the satellite and the rate of 
rotation of the earth around its axis are in a simple ratio. 

This problem has been examined from different points of view. Morando (1963) 
considered the particular, important ease of a 24-hour satellite; Allan (1967, 1973) 
worked out perturbations due to resonance in a more general case and to the first 
order, by using Lagrange’s equations. On the other hand, Garfinkel (1974) considered 
the abstract, mathematical problem of ’’ideal resonance." We have attempted to derive 
expressions for the perturbations in the motion of an artificial satellite due to the most 
general tesseral-harmonics resonance. The calculations are carried out to order 3/2 
in the small parameter of the expansion of the perturbing function, and the perturbation 
method used is that of Hori, which is canonically invariant and avoids mixing of old and 
new variables, as opposed to the generally used method of Von Zeipel. 


This work was supported in part by Grant NGR 09-015-002 from the National Aeronautics 
and Space Administration. 
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2. CHOICE OF THE DISTURBING FUNCTION AND EQUATIONS OF MOTION 

The gravitational potential of the earth at exterior points can be expressed as 
follows: 


U(r,6,X) = ^ 


« \ r 

n-2 


n 


oo n 


vn 


1 - \ ' J_. ( It) P n (cos 0) + EE ^n,m(r) ^n, ml COS ^ COS \i, 


m 


n=2 m=l 


( 1 ) 


where (r, 9, X) are spherical polar coordinates relative to the center of mass of the earth, 
the axis of rotation being the pole of coordinates; \i is the gravitational constant G times 


the mass of the earth; R is the mean equatorial radius of the earth; P r (z) is the n 
Legendre polynomial; P^ m (z) are associated Legendre functions: 


th 




/0 d n P (z) 

= n-9. 2\ m/2 2— 


dz 


n 


and J and J _ are dimensionless coefficients related to the normalized coefficients 
n n, m 

C , S— by: 
n, m’ n,m 


J n,m = J n,m V 2(2n + (n-m)J /(n + m)! , m*0 , 


- J n= C n,0' V^ C n,0 > 


( 2 ) 


J cos m(\ 1 = C ^ cos mX + S„ sin mX . 

TTl ' ti tyi 1 n tYi n rv» 


n, m' n, m 


It can be assumed that the largest perturbing forces are due to the leading zonal 
harmonic, which contains J ^ and corresponds to the oblateness of the earth. The 
remaining zonal harmonics have much smaller effects and are not considered here. 

The tesseral harmonics, containing , are the longitude-dependent terms and 
usually give rise only to short -period effects because they all contain the mean anomaly 
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and the sidereal angle in their arguments. In the case of resonance, however, a cer- 
tain number of these terms can produce effects of large amplitude and very long period. 
It is therefore necessary to select appropriately those tesseral harmonics to be 
included in our disturbing function. 


The general tesseral harmonic can be developed in terms of the osculating elliptic 
elements (a, e, I, £2, w, M), referred to the equator of the earth, in the form (Kaula, 
1966): 


n+1 

'n, m = (“Xa ) F n,m,p< I )2j G n,P,q (e) C ° S L (n " 2p) “ + <n “ 2p+q)M 

(p=0 q 

+ m ( n-vt-X n i ’ 


(3) 


where v is the angular velocity of the rotation of the earth, t is the time, and F (I) 

n, in ^ p 

and G n p ^(e) are, respectively, the inclination and eccentricity functions as defined 
by Kaula (1966). The functions p ^(e) are of order |q| in eccentricity. 


Resonance occurs when a pair (a, p) of mutually prime integers exists such that 
the satellite performs p nodal periods while the earth rotates a times relative to the 
precessing satellite's orbit plane. This can be expressed by 


a(co + M) = p(v- £2) 


(4) 


where to, M, and £2 are the rates of change with time of to, M, and £2, respectively. 
The corresponding slowly varying arguments are of the form 

k,$ + const , 

1 a,p > 


where 


4 = a(to + M) + P(£2- vt) 

a , p 

is called the resonant variable and ^ = 1, 2, 3, . . . . 
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Considering that d> is a small quantity and that the general argument in a tesseral 
harmonic is 

4> = (n - 2p) w + (n - 2p+ q)M + m{f2 - vt - X n m ) , 

we shall select the tesseral harmonics by keeping those containing arguments such 
that 


n - 2p + q = kjd 


m = kjP , 


k^ = 1, 2, 3, . . . 


(5) 


Since n s m and since lower order tesseral harmonics are bound to have larger 
effects because of the factor (R/r) n , where R/r < 1, it can be assumed that it is 
enough to consider only the cases k^ - 1, 2, and 3. 


If kj = 1, then m = (3, and the tesseral harmonics to be considered are 


V, 


P+k 0 ,(3 » 


kfj - Oj 1, 2 } . • . 


In a general manner, the tesseral harmonics to be kept are 


where kj = 1, 2, 3, . . . and for each k^ the index k Q takes values 0, 1, 2, . . . . 
2. 1 Resonant Part of a Tesseral Harmonic 


Once we have decided which tesseral harmonics to keep for our disturbing func- ’ 
tion, we then need to extract for each its most important "resonant" part. 
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If V is a selected tes serai harmonic, with n= k.B + k A and m = k.B, we look 
n,m * l r 0 l r ’ 

for terms containing "resonant" arguments by solving equation (5); that is, since we 
already know m = k^p, we have 

n - 2p + q = k^a , 

or, in this case, 

k^p-a) +k Q = 2p - q . 

As the condition 0 ^ p £ k^p + k^ must be satisfied, the following pairs (p, q) are 
solutions: 
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V = R + V' 
n, m n, m n, m 


} 


where V' is the residual, the effect of which is much smaller than that of R 

n, m 7 n, m 

It will be convenient to express the sum of the resonant parts of all the selected 
tes serai harmonics in the following manner: 


3 

2Xm = E Dfrj) c°s 2 0 i(ki) , 

k, =1 


where 


D(k.) cos 


Vi) - 


k o< k i> 




» 


the sum over k^k^) meaning that we have taken into account all the values of k^ when k^ 
has a given value. In order to obtain D(k^) and 0^{k^), we can write 


pR n J 




n, m 


n=k 1 p-t-k Q , m^kjp , 


A(k 1 ,k 0 , X )=F n)m?x (I)G njX>qx (e) , 


(8) 


x=n 


D(kp exp 2i0^{k^) = E S(k-^ ? kp) ^ ] A(kj> x) exp i(k^^ ^ - qui k^pX^ * 

k Q {ki) 


x=0 


where 


q x = k i (a -p) - k Q + 2x . 

If we write 

Ak p k o " klP V+VV 


and 
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then 


29 1 (k 1 ) = k^ p - k 2 (a - p) w + ip(. k a ) , 


D(k x ) exp ^(kj) 


n 

y s < k i» k 0 > ^2 a ^i» k o’ x) exp * [< k o - 2x ) w - A k k ] * 

k 0 (k x ) x=0 P 0 


Thus, 


x=k 1 P+k Q 


B{k 1 ) = D(k 2 ) cos ^(k A ) ^ S(k 1 ,k Q ) A(k p k 0 ,x) cos [(k 0 -2x)u - A fe fe J 

k 0 x=0 P 0 


(9a) 


x=k 1 P + k Q 


C(k 2 ) = D(k x ) sin ^(kj) = ^ ^ S(k p k 0 ) A<k p k 0 ,x) sin [(k 0 ~2x)to - fe 1 , 

k Q x=0 P ° J 


D(k x ) 2 = B^) 2 + C(k : ) 2 , 


(9b) 


(9c) 


in which D(k^) and $(k^) are functions of a, e, I, and £2. Appendix A gives expressions 
for the first and second derivatives of D (k ^ ) and i/'(k^) with respect to the canonical set 
of Delaunay variables. 


Finally, the potential in which the satellite moves is taken to be 


V = *± 
r 


J 2? 
J 2 3 
r 


P 2 (cose) + 2> njm+ £ 


v; 

n, m 


(10) 


the sums being taken over all the tesseral harmonics selected. 


2 . 2 The Equations of Motion 


If we consider the canonical set of Delaunay variables, 
L D = ^f :5: > G d = |/pa(l - e 2 ) , 


= M , 




H d = cos i , 

h D = , 


( 11 ) 
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then the equations of motion of the satellite are 


dL D 

8F d 

dG„ 

D _ 

8F d 

®D = 

9Fd 

dt 

“d ’ 

dt 

S «D ’ 

dt 


di_ 

D _ 

9F D 

dg D _ 

8F D 

dh^ 

D _ 

8F d 

dt 

" 8L d ' 

dt 

8G d ’ 

dt 

- 9 r d 


where is the Hamiltonian of the problem: 


F D ^72 + F 1 + F 2 » 

zlj D 


where 


R 


F 1 = - J 2 p 2 ( cos » 


( 12 ) 


(13a) 


(13b) 


which can be expressed in Delaunay variables, and 

F 2‘IXm • < ! 

Here, the Hamiltonian F^ depends explicitly on time. To avoid this, we can perform 
a canonical transformation so that the new angular variables will be 

l = £l> > g = §D ’ h=h D -vt . 

To find the new momenta and the new Hamiltonian, we have to solve 

L df + G dg + H dh - L D df - Gp dg - H D dh - (F D - F) dt » dV , 

where dV is the differential of a function. We keep the solution corresponding to 
dV = 0 as follows: 
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> 


L = L 


D 


G= G. 


D 


H= H. 


D 


F = F D + vH 


The equations of motion then become 


dL _ 3F 

dt ai ’ 


d£ _ 9F 

dt “ 3L » 


dG _ dF 
dt 3g ’ 


dg_ _9F 
dt " " 9G * 


dH _ 3F 
dt 15 5 


dh 9F 
dt “ 0H * 


with F = Fq + F^ + F^: 



(14) 


(15) 


r 2=E V n,m > 

in the above, u = g + f , and f is the true anomaly. 

The Hamiltonian is now expanded in terms of powers of the small parameter 


F x = 0(J 2 ) , 

*2 = °< J 2 > ‘ 
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3. HORI’S PERTURBATION METHOD BY LIE SERIES 


In order to remove the variables t and h from the Hamiltonian, we shall apply 
Hori's method by the use of Lie series (Hori, 1966; Aksnes, 1970). 


The greatest advantage of this method is the fact that the Poisson brackets that 
appear are canonically invariant, which means that each bracket can be calculated in 
terms of the most convenient set of canonical variables. Also, this method avoids 
mixing old and new variables. 


Let us consider a system of canonical equations 


_ 8F dy ,i ... 9F 

dt 3y. ’ dt 8x. ’ 

1 


where the Hamiltonian F is developed in powers of a small parameter c, 


(16) 


OO 


*-X> k , 

k=0 



and let us also consider a Lie transformation (Hori, 1966) 




x. 

1 


+ 9§_ + i j9S_ 
j ay 2 ) 9yf ’ 



, 9S_ _ 1 ( as 

y y j “ 8x( 2 \dx:’ 

J 1 3 



s| , s + . . . 
sj, s|+ ... 


7 


7 
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where { , } are Poisson brackets and S is a function of (x',y') of order 1 in a 


oo 



K = 1 


Then, if f{x, y) is a function of (x., y.), j = 1, n, it can be expressed in terms of the new 
variables (xf, y f) as follows: 


f(x, y) = f(x', y') + {f, S} + | {{f, S}, S} + | {{{f, S>, S}, S> + . . . . 


The equations of motion then become 





9F' 

9x' 

] 


where F' is the new Hamiltonian: 


(17) 




If we assume that F does not depend explicitly on time t, we can write the energy 
integral 


F(x,y) = F'(x',y') . 

Expanding in terms of e and collecting terms of the same order of magnitude, we 
obtain 


+ (*V ^1^ ^1 7 (^®) 

{ F 0- S 2> + ^ F l + r i> S l> + F 2 =F 2 ' 
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In order to average these equations, we introduce a pseudo time t': 

dx{ 9F 0 dy^ 9F Q 

dt/ &yT » dt' dxT ' 

J ) J 

If this system has a solution, then 

dS k 

< F o- s k> = -3r ’ k21 ’ 

and if A(t') is a periodic function of t', with period T, then 
T 

A g = I | A(t') dt' ( 19a ) 

0 

is the secular part of A and 

A = A - A„ (19b) 

p s 

is the periodic part of A. 

We can remove t' from F' by applying the averaging technique [eqs. (19)] to 
equations (18) (Hori, 1966; Aksnes, 1970). Then and F£ are uniquely determined 
by the following set of equations: 


S i = / F i P dt \ ■ 

F 2= F 2 S + I< F l + F i> S l>s - 

S 2= / F 2p + f{ F l + F i’ S l>p > 


( 20 ) 
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In the present case, 


F 2=Z D,k l )00s2e i (k l )+ E V n,m 


and the pseudo time t' is defined by 


di ' _ 
dt' 


dh' _ 
dt' 


!!_« jl. 

81/ ,,3 


f 



y 


so that 


dt' = 




( 21 ) 


In the process of calculating S^, we have to evaluate the following integral: 



) cos 20^(k^) dt' 


D(k 1 ) sin 20 '(kj) 
26 


> 


where 


6 = 


2 , 

3 ) a ~ V P| 

L>i /' f 


In the case of resonance of order (a, (3), 


n 


, JL-^ 

L' 3 Q 


} 


where n is the mean motion of the satellite; thus, 6 is a small quantity appearing in the 
denominator, and the preceding method cannot be used. 


14 



To avoid this situation, we shall introduce a procedure that is commonly used in 
different resonance problems (Garfinkel, 1973; Hori, 1960). The method consists 
of developing the Hamiltonian and the determining function S in powers of the square 
root of the small parameter: 

S = S l/2 + S 1 + S 3/2 + S 2 + * ' * ‘ 

If F' is the new Hamiltonian, then 

F'=F' + f; + F' /2 + F^ + ... 

and equations (18) now become 

r = F + {F 0 , S i / 2 } + {F 0 , S 1> + { F 0 ,S 3/2 } + { F i,Si> + {Fi, S i /2 } + ... 

+ 2 ^ F 0’ S l/2^ S l/2^ + 2^ F 0’ S 1/2^ ,S 1^ + 2 ^ F 0> S l^’ S l/2^ 

+ g {{(F 0 , S i/2^ S i/2^’ S l/2^ + 2 ^ F 0 J S l^ S l^ + 2 ^ F l> S l/2^ S l/2^ + " ' 

( 22 ) 


Equation (22) immediately yields to order 0: 

F^(L',G',H')=F 0 (L',G',H') . 

We shall group the different terms in equation (22) according to their order, but first 
we need to develop some of the Poisson brackets in the following manner: 

First consider 


o i _ P 2 9S l/2 . 9S l/2 

{ F 0> S l/2 } " ‘ ^3 a ~ + v W~ 


As Sjy 2 arises strictly from the existence of a resonance and the resonant terms are 
functions of the resonant parameter ^ = a (I '+ g') + (3h', we look for a solution 


S l/2 = S i/2< L '» G '> H '> %’> *[, p> such that 


15 



*1/2 d 8S l/2 

ai' (B ah' 

Thus, 


r F g x = ( JL. + }&\ 9S l/2 = aS l/2 
t F 0’ S l/2' L ,3 a j a' y dl' ’ 


2 3 

where y = <v|3/a)-(H /L' )is a small parameter that can be assumed to be of order 
1/2, so that {Fg, Si/g} is of order 1 in equation (22). 


Let us write 


S 1 =SJ 1 + SR 1 , 

where SJ^ arises only from the effect of J g , which can be obtained by solving 

F 'l - < F o> ®1> + F 1 =► j w,"- J l lp dt' - < 23 > 

and SR comes from the existence of the resonant terms that combine their influence 
with that of We shall look for a solution of the form 

SR X = SR 1 (L', G', H', g', p) 

such that 

9SR 

{ F 0 , SRj} = y gjp 

is of order 3/2. 
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We can then write 


Y = Y(L') , 


1 8y / 8S l/2\ 2 1 ( 9S l/2 

2 91/ \9i' / 2 Y )9i' ’ b l/2 


of order 1 


of order 3/2 


1 1 1 ( 9SR 1 

2 ^ F 0» S l^’ S l/2^ “ 2^ F 0 j ^1^ S l/2^ + 2 j Y aT ,S l/2 


2 {{F^SJ^S^} + 2 0L / d£ , 


1 ay 9SR i 8 S i/2 


+ 1 ( 9SR 1 s l 

+ 2 y | 9F“ ,S 1/2| 5 


2 {“t F o» ^1^’ ®l/2^ 


1 9y 9SR 1 9S l/2 

2 81/ at' at' 


are of order 3/2 and 


1 j^l 

2 Y ) at ' ’ s i/2 


is of order 2; 


jrr, F s > s 1 - 1 £* ( 0 ^/ 2 ) , 1 »V 

6 6 flT ,2 \a£ / / 2 ei/Var /) at' ’ 


6 y ))d£' > 0 l/2(’ 1/2 » 
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where 



is of order 3/2 and 



is of order 2. In the expressions above, 

o 9 9 

3y 3|i gy 12 jj“ 

^'l/ 4 ’ 8L' 2 "l' 5 ' 

We write Fg in the following manner: 

F 2 = 2 Z D(k l> si " 2 e < k l> - Z D(k l ) + Z V n, m < 

k l k l 

where 


2 ®l^l) > 

and D(kp is a function of (L', G', H', g'); we assume here that g' is reduced to its con- 
stant part gQ. The secular part of g' should be taken into account when the long-period 
perturbations are removed. 

We shall assume for the time being that our Hamiltonian contains terms for only 
one particular value of k,, and thus we can drop the k, index for convenience. 

We now proceed to splitting equation (22) by grouping terms of the same order. 

A. We shall allow the first approximation of the determining function, to 
contain the resonant part of the Hamiltonian, which means that it will contain secular 
terms. We then have, to order 1, 
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* F 0’ S l/2 } + I II 7 (aF^) + sin 2 0-0 , 


or 


Vfl/2) + 3 hL (fW 

\B£' / 2 l ,4 \aP / 


+ 2D sin^ 0=0 


This is called the "resonance equation. " 


If 


A - 2y^ 4 

a 2 y 
3|i ak^ 


k 2_ 12P^ 

t ,4 2 » 

L V 


( 27 ) 


(28) 


and we assume that = S^fS)* then by solving the quadratic equation (27), we get 


as 

ae 


^ = -A ± A y 1 - k 2 sin 2 ~0 , 


(29a) 


or 


as 

80 


1/2 _ 


= A (A - 1) 


(29b) 


if 

/ 2 2 

A - ± y 1 - k sin 0 . (29c) 

We discuss this in detail in Section 4. 
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B. To order 3/2, we have the following: 

aSRj ( 9S 
F 3/2 = Y ~dT~ + 2)m 


or 


1/2 „ ) 96 1/2 ^1 

7 » 1/2 j 8 L' a*' a F~ 


9S 0SJ 

+ 2 {{ F o> 81 J ’ S l/2^ + 2 1 l' aP dF + ^ F l> S l/2^ 

1 9 2 y / 9S 1/2 N | 3 1 &y 88 1/2 j 8S l/2 

+ 6 „ ,2 \9i' / 2 9L' 8 £' ) 9£' ’ 1 


9L 




F 3/2 yA 9^' " 2 )9f 


1/2 




1 1 rh 1 /2 1 

^ F l ,S l/2^ + 2 ^ F 0 ,SJ l^ ,S l/2^ + 2 9L' 9P 9i' 


(30) 


We notice that in equation (30), some terms contain only and are therefore due 
only to the resonant tesseral harmonics, while others express the interaction between 
the resonant terms and J^. Let us therefore write 


SR — Z ^ + Z 2 , 


where Z^ contains only those terms that are functions of S jy 2 > Z 2 contains all 
the others. If we write = X + Y, we then obtain two equations: 


X= yA 


+ s 1 1 9 2 V ( ^lM 

9 r 2 W ,s i/2j 6 / 


(31) 


and 


Y - s i/2> - 1 « F o> S V< s i/ 2 } + lit' sr 1 w 1 * 


9S, 9SJ i yA9Z f 

W 


(32) 
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C. Finally, to order 2, we have 


F 2 = F 2s = " D ( L '’ G '’ H '’ g,) > 


( 33 ) 


and 


0 " V n, m + ^ F 0> S 3/2^ + ^ F 0> S 2-^ + ^ F 0’ S 2^ 

4v|f?-. S l/2( + I {{F 0’ SJ l } > S, l } 


7; v||g/'' /2 > S l/2j'> S l/2.| + t F li S ll + 2 fl F l> S l/2l'> S l/2‘ 


+ ' 



L jh. (^/l\ 1 9S V 2 

4: |ai/\ai' / pF ’ 


as. 


i_ ay_ )) 00 1/2 


'as 


4 i ai/)} ar ,s i/2( ,s i/2( \ar 


1/2 


S l/2('» S l/2j 
(34) 


Equation (34) defines Sg, SJ,, and S^yg. The l erm ®2 ar ^ ses f rom the existence of 
J 2 only: 


{F Q ,S 2 } +|{{ f 0 » sj i}» sj i} +{ f p SJi} = 0 5 


(35) 


Sg, which takes care of the nonresonant tesseral part, is obtained by solving 


or 


{F n ,s;} + v' = 0 
1 O’ 2 J / > n, m 

S'=Y V' *' ; 

2 J n, m 


(36a) 

(36b) 
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and S 
of 6, 


3 ^ 2 is due to the resonant tesseral harmonics, 
such that 


which we assume to be a function 


^ F 0 ,S 3/2^ 


as 




3/2 


and 


. ^3/2 i !«]„ 

0 v a T + 2 Mat' ,s i/2( 2 8L' \a«' / at' 




+ 6 


1 ))8y 9S l/2 q l s ll jf 9Sl / 2 s > S 

6 nai? a a ,s i/2( ,s i( 6 n ai/ ar ,b i( ,s i/2l 


+ . . • 


(37) 


If the Hamiltonian contains more than one value of (usually kj = 1 is enough), 
we can write 




1/2 ^1/2, k, 

k, 1 


= ^J Zl ’ k i 

-?V. 


(38) 


We first obtain by solving equation (27), replacing by ^ for each 
value of kp and adding the different S l/2.k, terms. This amounts to neglecting mixed 
terms of the form ’ 1 
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^1/2, k x 9S l/2,k 2 

aF x W 


in this theory. 

We then solve equations (31) and (32) in the same way, again for each value of 
kj, and add the perturbations obtained; the neglected mixed terms are reasonably 
small. 


If (L,G, H, St, g,h) is the set of modified Delaunay variables as defined in Section 2, 
and (I/, G', H',i g', h') are the new variables, then after applying Hori's transforma- 
tion, the solution of the equations of motion (14) becomes 


L = 1/ + 


G = G ' + JiP + i {If 75 s } + • • • * 
H=H ' + w + i {w» s } + --' > 

1 = v ~ S? ~k ■ ••• » 

, as l (9S J _ 

S g 9G' ” 2 \9G° S f * ‘ ’ > 

v _ v, 9S_ 1 fas_ i _ 

_ 9H' ” 2 \9H" J » 


(39) 


where S = S(L / , G', H',i ', g', h') is the generating function 


S = ■*" SJ + SR + . . 


Then, if 




is the new Hamiltonian, we get 


23 



dl/ _ 3F > 

at a?' 

dG' _ 3F' 
dt £g' 


dH' 

dt 

dl ' 
dt 


d£l 

dt 


dh / 

dt 


®1=0 
a/ > 


2 9F' 


9D(k 1 ) 


§IL = jL _i + V _X = r 

" 31/ , ,3 31/ Z-< 81/ 1 1 ’ 


3F' 9F 1 9D(k P 

“ 3G' “ ” 3G' + Z-f 3G' S 1 * 

k l 

9F' 0F i v 

3H' " 3H' " v Z^ 3H' 1 ’ 

k. 


where I p g^, and h^ are constants such that 

i' = i't + i^ , 
g' = g\t + , 

h' = h't + h' Q , 

with l g^, and h^ also constants. 

In Section 4, we will give precise expressions of the perturbations up to order 
3/2. The perturbations due to J 2 only are not considered in this paper, but can be 
found in a classical theory of artificial satellites, for instance Brouwer's. However, 
the interactions between J_ and the resonant tesseral harmonics are expressed with 
the generating function Z^. 
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4. DETAILED EXPRESSIONS OF THE PERTURBATIONS 


In this section, we drop the primes in the variables (L', G', H', i g', h'), since 
no confusion exists here with the original Delaunay variables, which do not appear. 

4. 1 Expression of /„ 

As we already have seen, 

9S,/2 

ae =A(A_1) ’ 

where 

A = e |/'l - k 2 sin 2 0 , e = ± 1 , 


so that 

9 

S^ 2 = -A + eA j |/l - k 2 sin 2 x dx . (40) 

0 

Let us discuss the value of t in equation (40). First, we consider the case when 

2 2 

k < 1. For any value of 0, 1 - k sin 9 > 0, and we are in the circulation case. In 
order to have continuity with the classical solution, valid away from resonance, it is 
necessary that Sjy 2 — 0 when k — 0. Since A — » when k — 0, then, necessarily, 
c = +1. Otherwise, S^y 2 -co when A — <x>. 

For the case when k = 1, 

A = € cos 0 , 

0 

J A(x) dx = e Sin 0 . 

0 

In order to have continuity when k — 1 and k < 1, we must have e = +1. 
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When k > 1, 


i ,2.2 

1 - k sin x & 0- 


- arcsin ^ < x < arcsin ^ 


y 


and we are in the libration case. 


We need to know the expressions of the elliptic integrals of the first and second 
kind — F(0,k) and E(0,k), respectively - for the cases k > 1 and tt/ 2 < 0. For these, 
we have the following formulas (Gradshteyn and Ryzhik, 1966): 


E(mir ± 0,k) = 2mE ± E{0,k) , 

F(mTr± 0, k) = 2mk ± F(0, k) , 


( 41 ) 


where E and K are the complete elliptic integrals of 

E=E(f,k) , 

F=F(f,k) . 


Then, if k^ = 1/k, k' 2 = 1 - k 2 , and 0^ is defined by sin 0^ = k sin 0, we have 

E(0,k) = k : [k 2 E(0',kp + k' 2 F(0',k 1 )j , 

F(0,k) = k 1 F(0',k 1 ) . 


(42) 


Returning to our problem, we see that in order to have continuity when k — 1 and 
k > 1, we can take € = +1 in this ease also, because, from equations (42), 


lim E(0,k) = lim E(0', k) = E(0, 1) = sin 0 . 

k— -1 k— 1 

In conclusion, whatever the value of k, we can write Sjyg l 11 the form 


= -A0 + AE(0,k) . 
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To simplify this, we can set 


I 2 = E(0,k) , 

I 0 =F(9,k) . 


Then the derivatives of Sjy 2 are readily obtained: 


®®l/2 

8 r^=ac 2 Y(A- 1 ) , 

3S 

ag 1 ^ = dc 2 Y ( A + c 2 y£l 4 ^2 ” V ’ 

^= P o 2 Y(A- 1 ) , 

(43) 

9S 

9l” = f <V 9) + 8c l Y (I 2" 0) + c 2 p l Y(A_1) + < c 2 Ya l " 4 c 1 y > » 


®®l/o 

8G - c 2 y p 2 {A - 1) + c 2 Y a g (I 2 - 1„) 


J 


“l/fc 

9H 


= Y c 2 P 3 ( A - 1 ) + c 2 Ya 3 {I 2” I 0 ) 




where 


c 


2 



> 


and 





_ 1 8D 
a l D 9L 



9 


a = I 8D 0^ 

l D 9G ’ p 2 SG 


(44) 


[eq. (44) cont. on next page] 
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1 9D 
a 3 D 0H 


dip 

p 3 SH * 


1 3D = 8j£ 

a 4 D8g ’ P 4 8g 

d = (3 + p 4 . 


( 44 ) 


To calculate the elliptic integrals, we use the following expansions (from Gradshteyn 
and Ryzhik): 


0 < 0 < ^ , k 2 <l, 

F ( e,k)= 2 ( _ i /2 ) (-k 2 )“ t 2m (0) , 

m=0 

(m 2 ) ’ 

m=0 

where 

V6) = e > 

t 2 (0) = ^ (0 - sin 0 cos 0) , 

, ... 2m - 1 . ... 1 . . 2m - 1 _ 

W 6 * = "55” V - 1)< 0 > ' 2^ cos 0 Sin 0 ’ 

(-l) n (a) 

/a\ _ ' 'n 

\n/ nl > 

in which 


(45) 


(46) 


(a) n = a(a+l) ••• (a + n- 1) , n=l,2,... , 

Wo Bl > 
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and 


it k 

K(k)=|+ 2 ^ — m 

m= 1 


2m 


ml m! 


OO 




m = l 


2m) 



( 47 ) 


4. 2 Expression of and Its Derivatives 

We first calculate the Poisson brackets that appear in equation (31). After some 
algebra, 


( 9S l/2 c l „ L ,,/l ,2 ; . D . sinJcosG 

|9f ' ,S l/2j aC 2 Y ^Ha -1 )*^ <*0 A 

+ 4aC x Y (A-l) 2 + c 2 x y (A-1)(a--±) 

, /T T x ,2 sine cos 0 , „ „ ,,2 „ sin 0 cos ©1 

+ < I 2" I 0 ) c 2 YXk ~~ A + 4,10 l Yk (I 2' 9)— A j * 

(48) 


where c^, y, k, and A are as defined earlier, and 


,3 


c, = 


1 O 2 » 

3|i a 


x = | (da 2 + aa 1 + pa 3 -^a) , 


d = P + P 4 • 


(49) 
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The equation that gives X and is then 


8Z . qc 2 y A 

x ‘y A sr+— 2— 


-i) + k 2 (i 0 -e) 


sin 0 cos 0 


+ 4aC lY (A-l) 2 + c 2 xv (A-l) (a -^) + CgXyk 2 (Ig-ty 


sin 0 cos 0 


1 0 Y 3 2 /A- 1\° , . , 2 _ si 

+ 3 2 a C 2 y \—) + 4ac l^ k ( I 2" e ) — 

0L 


sin 0 cos 0 


( 50 ) 


Since (A- 1) is periodic in 0, we do not find any constant terms: 


X = 0 . 


We assumed earlier that SR^ was a function of 0, and we will now assume that Z is 
a function of 0 only. Therefore, 

9Z 1 _ 92 1 80 = a 

90 0f' 2 00 ’ 

and the following equation gives Z^t 



= ~ C 9^ 


(A - i)(i - 1) + k 2 (i 0 - 8) + 4aCjV (A _ 


+ c 2 x Y (A-l)(A-i) + c 2 YXk 2 (I 2 -I 0 ) 


sin 0 cos 0 


,, „ u 2 n sin 0 cos 0 4 „ ( K 2 0Al0 1\ 

+ 4aC lY k (I 2 -0) ^ 3 aC lY " 3A + 3 - - J 


(51) 


In order to integrate this equation and obtain Zp the following expressions are useful: 


30 



/ 

/ 

/ 


x sin x cos x dx _ _ xA(x) + 2 
|/l - k 2 sin 2 x 


I,(x) 


,2 ' ,2 * 

k k 


I 2 (x) sin x cos x dx A(x) I 2 j 
|/ 1 - k 2 sin 2 x 


k 2 + ]f \ x ■ T x * T 8in 2x ' ’ 


Ip(x) sin x cos x dx Iq A + 
^1-k 2 sin 2 x k 2 

Thus, by setting 


(52) 


q, = ac,Y 


and 


q 2 = c 2 yx , 


we obtain 


Zi=c 2 Y < I 2 (2 + q 2 + 8q 1 ) - 6 


s+ ? q i“T(f q i + 2q 2 ) 


sin 2e(“q 1 + 2q g ) + AI 2 (q 2 + 4q 1 ) - A0 (l + 4q 1 ) 




(53) 


Then the derivatives of with respect to the variables L, G, H, l , g, and h are 


9Z 

ax 


i = c 2 v| w [ 2 + «2 +8c ,1 + A <% + 4 ' 1 1>] 


90 

9X 

91, 


s + f Ol +2l l2) + A < 1 + 4q l> 


+ 8X t 1 q 2~3 q l + A(1 q 2*] 

+ fr [V^ + 4q i> - e <‘ + 4q i> + V 1 - q 2>] 

/20 „ \/k 2 90 

" V 3 q l + 2 q 2/ \ 2 9X 


1 r)k 

9X COS 26 + 4 9X Sln 20 


fclo / 2 k 2 

+ ^T I 2 + k 0 -T sin20 + AI 2“ I o- AI o; 


(54) 


[eq. (54) cant, on next page] 
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for X = 


with 


9 q, 


+ e -c 0 v 


2 Y 31 / 


20 k 4 

zu -- sin 20 + 4 AI 2 - 4 A0 - 1 1 0 


\ L_ 

+ ( 4 °lV + i-) |l 2 <2 + q 2 + 8 qi ) - 8 [ 3 + f q, (f q, + 2 qj 

- s i n 20 (t q l + 2 q 2) + AI 2 (q 2 + 4 q l> 

- A9 (1 + 4 qj) + I 0 (l - q 2 - | qj) + AI C (1 - q 2 |j (54) 
L, G, H, i , g, and h, where 


e = 1 
c = 0 


9A a x 


if 

if 


X= 1/ , 
X*L' ; 

2 


ax 


ax 


( A 1 \ k sin 0 cos 0 30 

( A "a) A 

8I 2 _ ae St 

A + 2 ^2 V 5 


ax ax 

ai 


a 


— — = — 4 — (X _ I \ 

ax a ax 2 V 


(55a) 

(55b) 


u 

t /q\ _ f dx _ 1 ( T ,2 sin 0 cos 0'\ 

-2 (0) - J ^3 - 77^ Kh - k A ) t 


*>2 

ax 


89_ 

ax 


(55c) 


„ ax , / 

3 2^ 8X + 6 ' 

% i \ 

4civ + -) x : 

(55d) 

30 P 1 

00 Q 


81/ 2 ’ 

df 2 ’ 


^ 80 _ P 2 
1 0G' 2 ’ 

39 _ P 4 

ag' 2 7 

(55e) 

80 P 3 

3H' ” 2 ’ 

89 £ 

ah' 2 • 
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X= l' or h' 



0 

4 2 

a l " L' a 

V 


if 

if X= L' 

if X = G' 

if X = H' 

if X = g' ; 


and 


ax 

ax 



aa, aa aa 2 ap 

ffir +p ffir + d ^r +i3 a3r 



(55f) 


(55g) 


4.3 Expression of Z „ and Its Derivatives 

From equation (32) and from the fact that {Fq, SJj} = - F^ , we can derive Z 9 as 
follows: 


1 9Z 2 
y = {f p S !/ 2 } - 1 {f i p > s !/2 } + yA ar 


1 8* 9Sl /2 8SJ 

2 9L 3i dt 


(56) 


Here we can neglect Fj , the periodic part of F j, as it is of order 1 in eccentricity 
and therefore very small for the satellites we are dealing with. We then have 

( ) 9 ^2 

Y =fi S ’ s i/2/ +vA ar • < 

Since there are no constant terms, 

Y=0 , 


and therefore. 


az r 


as 


y a 


a i 


■= - lb. 


1/2 


as 


+ b 


i/2 


as 


2 ag 


+ b 


1 / 2 ’ 


3 ah 


(58) 


where 
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9F ls 

4 2 / 

3J 2 M R / 

■ -i-f 3 

4 4 

H 2 

3L 

g 3 l 4 V 

G 2 

9F ls 

3J 2 ^ 4 R 2 f 

1+5 

h_ 2N 

9G 

1 

It 4 

CO 

O. 

4 4 

G 2 / 

0F, 

Is 

4 9 

3 J 2 fi R^H 



3H 

+ 2 l 3 g 5 

? 



and, since is assumed to be a function of 0, 

3^2 2c„ /A - 1\ ( 

TXT’-— |< V + V + b 3 d > (V) + b 2 a 4 V 


h - [ 0 


If we write 


B 1 - - — < b l“ + b 2 d + b 3« 

2 Co 

B 2*- — b 2 a 4 > 
we can obtain from equation (59): 

e , , , 


f I 2 (x) - I 0 (x) 

z 2 = B l( e-i 0 ) + B 2 J 

0 


Now we can write the integrals (Gradshteyn and Ryzhik, 1966) 

e 


J - 

J Ai 


(X) I 


— — dx= — 
A(x) 2 5 


fyx) _i„e 

A(x) ^ 2K log 


®(I 0 ) 
0 ( 0 ) ’ 
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where E and K are complete elliptic integrals of the first and second kind, 
and 0 is the following theta function: 

“ 2 

0 (u) = 1 + 2 2_] qm cos < 2mv) » 
m=l 


in which 


/-irk^ 
q = exp J 





and 


v = 


TTU 

2K 


The derivatives of Z^, then, are 

az 


ax 


2 _ 1 89. fo\ . . t i 

~ B, ( - rj-v ) + aY ^0 


6B 

'l V ax " ax/ ' ax 

®2 f VVg 
ax J a(x) 

0 


ae (h^o 


+ i ~ , — — dx + B 

£a (7 A 

.2 _ , ,2 


2 ax V A / + B 2 5 


where X= L', G', H', l', g',h' and, if k' = 1-k , 


S 

11 

1 

a 

X 

r*k’* t 

0 . + ~ - I 

2k' 2 

2 2 V 

L- 

8B 1 

2c 2 / 

8b 1 ab 2 0 

ax 

— V 

ax a+ ax d + £ 

8B 2 


/dh 2 + . . 9a 4 

ax 

— 1 

vax a 4 d 2 ax 


2 2 ' 
k sin 8 cos 9 , k sin 

A + S~ 


9P 


ab. 


where 


respectively, 


(63a) 


(63b) 


(64) 


9 


(65a) 
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€ = 0 if X * 1/ 


= 1 if X = L' ; 


II 

0 if X = i ' or 

9p 4 

7x7* 

9p 4 ep 4 0 p 4 

8G 7 ’ W* 311(1 ag 7- 

ab : 
ax ' 

ab 2 ab 3 . .. 

ax ax 0 lf 

ab A 

4b i 

81/ ~ 

L’ 

ab„ 

3 h • 

aT7~ 

_ L' b 2 ’ 

ab 3 

ai7 

3 , 

' L' b 3 5 

8b i 

3b l 9 J 2 I j4r21 

"8G 7 = 

G ' + 2 L' 4 G' 

ab n 

2 

4b l 15 

0G 7- 

g' + 2 l ,3 g ; 

ab 3 

5b 3 

3G' 

" ~G r ; 

8b x 

9 J 2 p 4 R 2 h' 

air" 

2 L/ 4 G' 5 

ab 2 

15 J 2 M 4r2 H'- 

BIT 

2 l ,3 g ,6 ’ 

ab 3 

t 4_2 

3 J 2^ R 

air 

2 L' 3 G' 5 


are given in Appendix A ; 
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We also need the second derivatives of Z g that appear in expressions such as 
{az 2 /ax, S^} for X= L',G', H',i',g',h', as they are used in the final expression of 
the perturbations. Their equations follow: 


where 


0 2 Z, 


9B, 


9X 9Y 9Y 



9 2 B, 


+ l*W< 0 - I o> + B i 


9 2 0 


6 2 I 


0 


9 2 B, 


9X 9Y 9X 9Y 
0 


9X 9Y 


f 1 2 X 0 J 0B 2 /}„ 9B 2 90 f 1 2 " 1 ! 

J A ^ + 9X t&W 9Y 9X \ A 
0 


+ B 
9B 


a 2 e (hjo) +-a 

2 9X 9Y V A / 2 9X 


(9I 2 /0Y) - (9I 0 /9Y) j 9A 


0 


9Y 


i^(X) + B 2 J* 


9Y 9X 


*2 I 0\ 1 


-72 

A 


( 66 ) 


90 _ 1 9 0 

9X 0Y 2 9X 9Y 


if X and Y are L', G', H', or g', 


[see eqs. (A-9) in Appendix A] 


9 2 9 


9X 9Y 


- 0 if (X or Y) is (i ' or h') ; 


9 2 B, 
9X 9Y 


2 c 


(z\ 


a 2 b. 


2 ( w W 1 “ "2 9b 2 9p 4 ^4 9b 2 

a~ \9X 9Y a + 9X 0Y d + 9X 9Y + W 8Y~ 


in which 


aV 


9 2 b, 




+ 

9X 9Y b 2 






+ T) 

a 

V§x a + 

6 = 0 

if 

X* 1/ , 

e = 1 

if 

x = 

= L' , 

TJ = 0 

if 

Y * L' , 

V- 1 

if 

Y = 

= L' ; 


\ 9B , 

W “ 17 0Y~ 


9b, 


9X 


2 ^4 9b 3 \ e 

d + wr b 2 + gx 7 72 B 1 
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a 2 b, A- 9 2 b 3 

ax 9Y = 9X 9Y = 9X 9Y = 0 lf < X or Y > = V or or 


8 b l _ 20 
9L' 2 ' L' 2 1 

8 b 2 _ 12 . 
9L' 2 L' 2 2 


9 2 b, 


12 


9L' 2 L' 2 3 


9 2 b 


9 2 b, 


4 8b l 


1 _ v “1 
3l/ 9G' 9G' 91/ ' T7" 0G 


0 2 b, 


9 2 b, 


4 8b l 


91/ 9H' 9H' 9L' " 1/ 9H' » 


9 2 b x 
9G' 2 G 


3 u 3 ab I 27J 2 m 4 H 2 H' 2 
,2 b l G 7- 9G 7 ” l ,4 ,7 


9 2 b, 


9 2 b, 


45 


J 2 P 


4 R 2 H' 


9G' 9H' 

9H' 9G' 2 

L' 4 G 

f5 

a 2 b x 

9 J 2M 4 R 2 



9H' 2 

2 L' 4 G' 5 




_ a2b 2 3 

9b 2 


91/ 9G' 

9G' 91/ V 

9G' 

J 

9 2 b 0 

2 

d \ 3 

9b 2 


91/ 9H' 

9H' 91/ “ ‘ 1/ 

9H 7 

> 

a 2 b 0 

2 

105 4 2 ^ R 

4 

9b 2 

9G' 2 

~ L' 3 G' 8 

G' 

9G' 


G* 2 "2 ’ 


(67) 


[eq. (67) cont. on next page] 
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b\ 

e 2 b 2 


15 J 2 ji 

8G' BH 1 

T BH' BG' 

G' IST 

L' 3 1 

a 2 b 2 _ 

4 2 
15 J 2^ R 



BH' 2 

2 l ,3 q ,5 

> 


8% 

e\ 

3 8b 3 


81/ 8G 

' 8G' BL' 

L' BG' 

> 


3 2 b 3 

3 ab 3 

• 

81/ BH 

T 8H' 3L' _ 

L' 9H' 

♦ 

9 2 b 3 _ 

30 h . 



8G' 2 

O' 2 3 ’ 



a 3 

_ d \ 

5 3b 3 



9G' BH' BH' 9G' " G' BH' 5 



8 % j ± 

ax 3Y ax BY 


/_L\ 1 BA dijj 

\2AS ~ OA 2 BY aX 
2A 


l^x 

+ 2 BY 


^-2 ” V 


(67) 



o O 

3k sin 9 cos 0 , k sin 9 cos 0 BA 
W A + a 2 BY 


- k 


2 c os 20 
4A 


8 i/> 

BY 


k'^ 


8 k 

BY 




k sin 0 cos 0 


b 2 b 2 2c 2 /a 2 b 2 ab 2 aa 4 ab 2 aa 4 \ 
ax aY = ' IT \ax by a 4 + W W + by ax / 



+ v 


8C 1 f^V +b !VL_S_b I 

” \8X" a 4 + b 2 BX / l ,2 B 2] • 
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2 

Note : In the expression 9 Z^/dX 9Y, we have neglected the term 



0 


which is a very small quantity. 


We also give the second derivatives of S.^, so that we can calculate Poisson 
brackets: 


| 9S l/2 g 
9X ’ d 1/2 f > 


, 9S 


1/2 


i ax 


’ Z U 


9 S l/2 z 
I0X ’ 2< 


• * • • 


If the first derivatives of Sjy 2 are in compact form, 


as 


1/2 


, 2 ae 
ax ^ c 2 y Vax ax 


ai 


+ .^-61(80^ + 1) , 


where 


then 


€ = 0 if X ± L' 


and 


e =1 if X=L' , 


A i/2 „ ( a \ s 2 e ' 

ax 8Y 2 y \ ax 9Y " ax 9Yy 


-&- i ) Kv + D 


/2 


( 8 c i Y + f ) (5 - i) + 86 A - 0 ) Y + ^ 


where 


and 


r7=0 if Y * 1/ 


a\ 


and 


tj = 1 if Y = L' , 


.1 \ tu A^A 0A Aft 

ax 9Y 2 0 Y ( 2 V 2 9 Y 2 SY 2 3Y 8Y J 
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5. REMOVAL OF LONG-PERIOD TERMS 


If we do not assume that g' is constant, as we did in Section 4, then, after 
eliminating t ' and h', we are left with the following system of equations: 


dL' _ 

8F 

df ' 

8F 

dt 

8i ' » 

dt 

“ w 


* 

dg^_ 

* 

l» 

11 

8F 

8F 

dt 

9g' » 

dt 

" 8G' 


* 


* 

dH' 

8F 

dh' 

8F 

dt 

8h' ’ 

dt 

~ 8H' 


where 


* * * * 

F = F 0 + F 1 + F 2 > 


F 0 = J=L 2 +vH ' ’ 

0 2L' 2 

F*=F*(L',G',H') , 

F * = - ^ °( k i> = F 2< L '’ G '> H '> g,) 


( 68 ) 


(69) 


We now proceed to remove g' from the Hamiltonian, by performing a new canonical 
transformation using Hori's method. 

** 

We introduce a determining function S' and a new Hamiltonian F such that, if 
L", G", H", l ", g", and h" are the new variables. 



and 


_ * - * ** 
{f c ,s'} + f 2 =f 2 


^ * 

Introducing a pseudotime t , defined by 

dg« 

u fr ~ i /t" G" H"^ 

** 9G y/ ^ ^ 9 n ) 

dt 


or 


v-1 


dt** = pdg" 


with p - \ - 


aF l 

dG 77 / 


** 


we can solve for and S' * 


** * 

^2 = F 2s ’ 
r * , * 
< F 0> S '}'- F 2P 


■/ 


* ** 

F 2p dt 


** 


Let us find the secular and periodic parts of F, 

l 

fixed value of kp the complete expressions of F 2 




By using expressions for only one 
and S' are then 


** 


F 


2 



k 


1 


** 

F 2 5 


s' = S^) • 

k. 


We can assume that is fixed and can drop the index k^ in the following derivation. 
We then have 

* 

F = -D 
2 > 

where 
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in which 


Aj - S(kp kg) A(kp k^, x) , 

*x'( k 0 - 2x ) g -\, k(( > 

and the summations are taken over values of x and [see eqs. (8) and (9)]. Hence, 


d2 =E a i - 2 E Vs - 4 E A i A t sin2 (V - 1 


i<] 


i<] 




sin [(<j).-<t).)/2]| 


To the first order. 


-4>j)/2]| 


D.fr A - 

1 ■ h (z>j 

Y'aA. Va.A. cos (*.-*.) 

D = y" A .-^-i + ^^ — -i— ai 

t 1 Z)a. Ea. 


Therefore, 


and 



(70a) 


(70b) 
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and, if X = L", G", or H", 
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The new Hamiltonian becomes 


** ** ** ** 

F = F 0 +F 1 tF 2 - 


** I ! 

F = r ^ 
0 2L" 2 


+ vH" , 


F 1 = F ls <L",G",H") , < 7 ' 

3|C 

F 2 = F 2s (L" f G", H") ; 

and the passage from the old variables (I/, G', H', i g', h') to the new ones is given 
by 


T ' - T " 4- M + 1 i g ,l 

L_L 9i" + 2 / *“ ’ 

c'- C" + + 1 |i!§l S '\ + 

G - G +9^77 + 2 W” S J + • * * » 


V = i 1 


as' 

1 

(as' 

~W> 

+ 2 

tah 77 ’ 

9S' 

1 

las' 

9L" 

“ 2 

\aL" } 

9S' 

1 

/9S' 

9G" 

2 

\aG"’ 

9S' 

1 

/as' 

9H" 

2 

\ 9H"’ 


The new system of equations is given by 



** 

di"_ 

** 

dL" _ 

9F 

aF 

dt 

dt" > 

dt 

“ 9L" 


** 


** 

dG" 

9F 

dg" _ 

9F 

dt 

ag" ’ 

dt 

9G' 7_ 


** 

dh" 

** 

dH" 

aF 

9F 

dt 

ah" » 

dt 

9H" 
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and, since F = F (L",G",H"), the solution is 


L" = L" = constant , 
G" = Gq = constant , 
H" = H" = constant , 
t"=t" + £", t , 

g"=g^ + g",t , 


’0 

1" 

0 


h"=h" + h",t , 


*o» 

h" i " 

V l \> 

<*" h" 
*»1» n l 


** 

2 

n _ 

3F 

± 

1 

" 8L" 

L" 3 


** 

8F 

tt — 

3F 

* Is 

1 

"aG 77- " 

9G" 


8F ls 0F 2s 


SF 


2s 


** 


8G" ? 

* 


aF _ 9F ls 0F 2 s 

h l ’W r ~ ~ v ~ ~M ir 


3H" 


( 75 ) 


(76) 


Returning to the initial set of modified Delaunay variables (L, G, II, H , g,h), we 
now have 

L = L' + §r (L', G', H', i g', h') + 1 {fr , s} + . • . , 

G = G' + (L' , G', H', l g', h') + |{J|r , s} + - . • , 

H “ H ' + fr <L', G', H', 1 g', h') + | {fy , s} + . . . , ^ 

i = i ■ ' - H? (L', G', H', V , g', V) - 1{||7 , s} - . . . , 

S = S' -|§ 7 <I/,G',H',i',g',h') -•|{|| 7 ,s}- ... , 

h = h' - |§7 (L', G', H', t g', h') - \ {||7 , s} - . . • , 

where S= S^ 2 + SJ + + Z g , which, with the help of equations (28), gives the com- 

plete solution. 
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APPENDIX A 


DERIVATIVES OF D(k x ) AND jp(k ) 

A. 1 First Derivatives 


We recall the following equations from Section 3: 

D 2 (k 1 )=B 2 (k 1 )+C 2 (k 1 ) , 

\ 

D(k : ) cos </»(k 1 ) = Bfkp , 

D(kj) sin ^(k^) = C(k 1 ) . 


(A-l) 


Then, if X is any of the Delaunay variables L, G, H, or g, 


D — "R 4. p QC 

ax _ B ex c ax 4 


BXfg, 


kiP+k 0 


ax 12 ^ 2 i ax A < k i> k o> x) ] cos [ (k o 2x,g A k , , k A ] 

U « . — A A U 


k 0 x=0 


8C 

ax 


k i^ k o 

12 12 ft [ S < k l' k 0> A < k l> k 0’ x) ] sin [<V 2x)g - Afcj.kJ > 


k Q x=0 


aA(k 1 ,k rt ,x) 3F (I) _ T 3G (e) „ 

IL 0 = 5j. m ;P 9L G , e \ + F m ibP.q 1 [ 9e_ 

ax ai ax U(e) + n,m, P (I) ae ax 5 


(A-2) 


asfl^) ^ o 

”2 (pk 1 +k Q + 1) 


if X * L , 


ax 


L S(k p k 0 ) if X = L 


A-l 



Ifx=g, 


l* k i +k 0 


9 B _ 
3 g 

-2 

2 


k o 

x =0 



pk : + k 

9 C 

3 g“ 

2 

2 


k o 

x =0 


— S(k|, k^) A(kj, ^qj *) {k 0 - 2 x) cos (kQ 2 x)g 


A, 


For the derivatives of i/>(k^), we have 
C(kj) = D<k 1 ) sin ^(kj) , 

so that 


3 C(k,) 


80 (k.) 


aDO^) 


OX » D(k 1 ) -jgg — cos ip( kp + -gg sin ip( kp 


and 


3 ^(k,) 


9 X B(k.) 


9 C(kj) x 9 D(k 1 ) 


9 X D(k 1 ) 9 X 

As we did in Section 2 , we can now write 


C(k.) 


a - 1 90 nr > 

a l D 1 (k 1 ) 9 L 'V * 


1 9 D 

3 D(k 1 ) 9H (K l' ’ 

p x = ■§!: < k l> » 

P 2 = II (k l> » 


a = -±- ® (k ) 

2 D(k 1 ) 9 G 1 1 ; > 


a 4 D(k 2 ) dg ’ 

P 3 = " 9 H 9 


_ M 


4 9 g 


(k,) 


i* k J 

(A - 3 ) 

' k] * 
1 * K 0 J 


(A -4 a) 


(A - 4 b) 


A -2 



A. 2 Second Derivatives 


Let us drop the index for simplicity: 


I 3D- 1 
d ax d 2 


B 


SB 

ax 


+ c 


ac 

ax 


a_ (i ao\ 
aY \p ax/ 



i 9D r ®, c 9C 
D 3Y A 8X ^ ax 


i / 3B aB ac ac r a 2 c a 2 B \ 

+ d 2 \sy ax aY ax ax aY ax aY / 


(A -5) 


If Y + g and X * g, then 


lrW=£ 

k o x 

2 2 

ax°aY 51 ^2 ^2 ax aY (SA) sin [ (k o " 2x)g _ A k p kj 

k A x 


(A -6a) 


(A -6b) 


If X # g and Y = g, we have 


and 


a 2 B _ V 9SA 

ax ag Z-j 2~t ax 


(k 0 -2x)sin[(k 0 -2x)g-A^ m ] 


a 2 c 
ax ag 


V V' a§A 

Z-j Z^t ax 


(k Q - 2x) cos [(k Q - 2x)g - A^ m ] . 


(A -7 a) 


(A -7b) 


A -3 



If X = g and Y = g, then 


and 



HE a(k o ' 2x)2 C0S K ■ 2x)g ■ *k lf kj 

k rt x 


(A -8a) 



EE 54 * 11 #' 2x)2 sin K - 2x » - k 0 ] 

k n x 


(A -8b) 


We find the following second derivative of Jp(k^) from equations (A -4): 


9B M 4. -r ^ 1 9D r 9_ (l BOX _ 9 2 C 

9Y dX a ax 3Y 9YD8X 9Y \D 9X j 9X 9Y 


(A -9 a) 


so that 


9 2 ^ 

_ 1 

9 2 C 

9B 9^ 1 3D 9C 

c ~ 

(l 

9X 9Y 

B 

9X 9Y 

" 9Y 9X " D 9X 9Y ~ 

9Y 

\D 9X J 


(A -9b) 


A -4 
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